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We combine the constraints on the pion quark structure available from perturbative QCD, non- 
perturbative QCD (nonlocal QCD sum rules and light cone sum rules) with the analysis of current 
data on (Q 2 ), including recent high-precision lattice calculations of the second moment of 

the pion’s distribution amplitude. We supplement these constraints with those extracted from the 
renormalon approach by means of the twist-four contributions to the pion distribution amplitude 
in order to further increase stability with respect to related theoretical uncertainties. We show 
which regions in the space of the first two non-trivial Gegenbauer coefficients 02 and 04 of all these 
constraints overlap, tagging this way the pion structure to the highest degree possible at present. 

PACS numbers: 12.38.Aw,12.38.Cy,12.38.Lg,13.40.Gp 


The pion distribution amplitude (DA), ip n , i plays 
a key role in hard-scattering QCD processes because 
it encapsulates the essential nonperturbative features 
of the pion’s internal structure in terms of the par- 
tons’ longitudinal momentum fractions Xi. Many studies 
11111111 have been performed in the lit¬ 
erature to determine this DA using the high-precision 
CLEO data |Tol| on the pion-photon transition, or to in¬ 
fer it from the calculation of the corresponding form fac¬ 
tor F VJ *ry(Q 2 * ). In particular, in | 8 j we have used light 
cone sum rules (LCSR) im with a spectral density ob¬ 
tained in the standard factorization scheme at the next- 
to-leading-order (NLO) of perturbative QCD. We exam¬ 
ined the theoretical uncertainties involved in the CLEO- 
data analysis in order to extract more reliably the first 
two non-trivial Gegenbauer coefficients 02 and 04 , which 
parameterize the deviation from the asymptotic expres¬ 
sion ip^ sy = 6xx (where x = 1 — x). A major outcome of 
this study js] was that the best agreement with the data 
is provided by a pion DA with a 2 ® st fit (/ig Y ) = +0.23 
and a 4 ® st fit (/ig Y ) = —0.22 at the scale /ig Y = 5.76 GeV 2 
(Schmedding and Yakovlev Q), where the twist-four con¬ 
tribution was represented by its asymptotic form. This 
best-fit solution (+) turns out to be very close to the 
model we derived before in [Q dubbed BMS model 
(x) —from QCD sum rules with nonlocal condensates 
(NLC QCD SR)—see Fig. Q Note that this model DA 
yields a prediction for F nj * 7 (Q 2 ) that is in good agree¬ 
ment 0, [3 not only with the CLEO data, but also with 
the older low-Q 2 CELLO (u| data. The profile of this 
pion DA is double-humped, but has its endpoints x = 0 
and i = l strongly suppressed. This suppression is di¬ 
rectly related to the vacuum quark virtuality A 2 , whose 


value 0.4 GeV 2 , extracted from the CLEO data in Q, 
coincides with that used in [ 13 ], the reason being that, in 
the latter, it guarantees best stability of the NLC QCD 
SRs employed. As one sees from Fig.0 the region of the 
02 and 04 values, admitted by the NLC QCD SRs (small 
slanted rectangle), is almost entirely enclosed by the lcr- 
error ellipse. The exclusion of the Chernyak-Zhitnitsky 
(CZ) pion DA (■) and the asymptotic one (♦), found 
before by Schmedding and Yakovlev Q at the 2cr level, 
was reinforced in Q. In the updated analysis Q, which 
accounts for the variation of the size of the twist-four un¬ 
certainties more precisely, these pion DAs were found to 
be outside the 4cr and 3cr ellipses, respectively. 

Recently, high-precision lattice measurements of the 
second moment (£ 2 ) 7r = fg(2x — l) 2 <p n (x) dx of the pion 
DA were reported by two different collaborations I17lll8j . 
Both groups extracted from their respective simulations, 
values of at the scale /i| Y , but with different error 
bars . 1 Unfortunately, the determination of the value 
of a 4 via the moment (^ 4 ) w appears to be a very diffi¬ 
cult task on the lattice because it involves four covari¬ 
ant derivatives E 2 - It is remarkable that these lattice 
results are in striking agreement with the previous esti¬ 
mates of 02 both from NLC QCD SRs [l2| and also from 
the CLEO-data analyses—based on LCSR— IBS , as 
illustrated in Fig. ^ where the lattice results are shown 
in the form of a vertical strip, containing the central 
value with associated errors. The (preliminary) results 
of the QCDSF/UKQCD collaboration da are shown in 
the left panels of Fig. Q The right panels of this figure 
show the latest estimates by Del Debbio et ah, reported 
in 0 . Previous determinations of the same team M 
favored the asymptotic DA, on the basis of much smaller 
values of the second moment (£ 2 ) 7 r. It seems that the 
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1 We have evolved at the two-loop level both lattice results to the 

scale /Xg Y in order to facilitate comparison with the CLEO-data 

analyses. 
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FIG. 1: The top figures show the LCSR-based CLEO-data analysis |j), 3 on F-n- 1 *-y{Q 2 ) in comparison with the lattice results of 
03 (left panel) and M (right panel). The displayed models are explained in Table [I] The bottom figures show the analogous 
results for the present data analysis which includes twist-four contributions by means of the renormalon approach. Here we 
use: dashed line: lcr; solid line: 2a; dash-dotted line: 3a); slanted shaded rectangle: confidence region of NLC QCD SRs m 
for A 2 = 0.4 GeV 2 . All results are evaluated at p,§ Y = 5.76 GeV 2 after NLO ERBL evolution. 


higher the precision of their simulation, the larger the 
result for (£ 2 ),r and, consequently, the more the asymp¬ 
totic pion DA is disfavored. Remarkably, the value of 
02 of both displayed lattice measurements (middle line 
of each strip) is very close to the CLEO best fit in (j. 3 
(+), whereas the model pion DA (Jt) 12] is either in¬ 
side the error margins (Del Debbio et ah), or just on the 
boundary (QCDSF/UKQCD collaboration). 

But the success of the CLEO-data analysis depends 
crucially on the progress in reducing theoretical un¬ 
certainties originating from uncalculated (perturbative) 
higher-order radiative corrections and (mainly nonper- 
turbative) higher-twist contributions. Those next-to- 
next-to-leading (NNLO) corrections proportional to the 
lowest /3-function coefficient in the standard scheme were 
recently calculated in f20] and addressed in Q. The 
higher-twist contributions encode information on power- 
suppressed corrections, related to higher quark-gluon 
Fock states and the quark transverse momentum. Just 
recently, Agaev 0] has attempted to improve the CLEO- 
data analysis by including twist-four contributions, spec¬ 
ified by the renormalon approach [22:, 23, 24} ■ The main 
aim was to estimate an upper bound for this contribution 
and extract new constraints on the Gegenbauer coeffi¬ 
cients 02 and 04 of the leading twist-two pion DA. Agaev 

I 


claims that the renormalon approach excludes negative 
values of 04 . This finding, remarkable if true, is in¬ 
compatible with the results derived from the analyses 
in @,1.7, 8, [13, discussed above, and, therefore, deserves 
closer scrutiny. 

Retracing Agaev’s calculation, we found that it is seri¬ 
ously flawed in several respects. Firstly, we revealed that 
the normalization of the twist-four pion DAs is wrong 
because it erroneously depends on the Gegenbauer co¬ 
efficient 04 instead of being a constant. This might be 
the root cause for excluding negative values of 04. Sec¬ 
ondly, the Gegenbauer coefficients evolved to the scale 
H g Y , given in his equation (4.2), turn out to be incon¬ 
sistently low. These errors aside, this sort of approach 
looks quite attractive to be followed and rectified. In 
an effort to carry the renormalon approach to its logical 
conclusion and compare its predictions with the previous 
ones, we take recourse in the present inve stig ation to the 
robust results of Braun-Gardi-Gottwald 24|, which gen¬ 
eralize the previous work of Andersen in |22l ], and embed 
their twist-four contributions to the pion DA into our 
data-processing framework @ @ ■ 

The starting point of our discussion is the twist ex¬ 
pansion of the matrix element that defines the pion DAs 

(e.g., 0 ): 


( 0 | |"d(z) 7 J/ 75 >V(z, 0 )m( 0)1 |tt + (p)) = if n dxe lxpz \p„ |V 2 ) 0 ;/ 4 ) + z 2 <p[ 4) (x; pi) 

L V! J 0 l L 


+ ( z v (pz) — PvZ 2 ') t4 4) ( x > Mf) } + °( z4 ) - 


( 1 ) 
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where the scale fip defines the factorization scale. Gauge 
invariance is ensured by the insertion of the light-like Wil¬ 
son line W(z, 0) = Pexp — ig f^dt z M A M (z — tz) . The 
twist-two contribution to the pion DA can be represented 
as an expansion in terms of Gegenbauer polynomials to 
read 


^ (2) 0) = a n^n{x) ; V’n(z) = 6x2 • C?J 2 ( 2 x - 1 ) . ( 2 ) 


n —0 


For the leading twist-two part, we use a two-parameter 
model with the coefficients 02 and 04 . This model re¬ 
ceives support from nonlocal QCD sum rules El (cf. Fig. 
0 . At the reference scale /rg Y , we have after NLO evo¬ 
lution af MS = +0.14, a 4 MS = —0.09. Note that though 
the NLO evolution generates higher Gegenbauer harmon¬ 
ics in the pion DA, the corresponding coefficients remain 
numerically negligible, so that the analysis can be based 
only on 02 and 04 @J. To estimate the twist-four contri¬ 
bution, we now use the renormalon approach of Ref. (3 ■ 
To this end, we express the corresponding DA, d>( 4 )(x), 
characteristic for the reaction 7*7 —> n° ;U! , with the 
aid of (x), in the form of a convolution 

$ ( 4 ) (x;/i 2 ) = 5' 2 (n 2 )K(x,y)®<p {2 \y) 

= 8 2 (p 2 ) f K(x,y)ip (2 \y)dy , (3) 

Jo 

I 


where the kernel K is determined in the present work to 
be 


K(x,y) 


|{% > x) 


XX 1, / X 

— + - In ( 1 - - 

IV V 


+ (a 


x,y 



(4) 


In Eq. 0 the coupling 8 2 (p 2 ) is defined by (see, for 
instance, m (n(,p)\g s dG afl j a u\0) = where 

G ap = (l/2)e aflp crG pa and G pa = G% A“/2. The nor¬ 
malization equation for the twist-four contribution, fol¬ 
lowing from Eqs. 0 and 0 , is 

r 1 g 

/ dxK(x,y)®tl) n {y) = -5 0n . (5) 

Jo y 

It is important to notice in this context that Eq. 0 ac¬ 
cumulates the content of the twist-four pion DAs, termed 
+ 2^1 4>||, 4/|| in 0 (Eq. (3.4) there)—see also 
Eq. (20) in 111 . Explicit expressions for these DAs can 
be found in l24l| (Eqs. (3.4) and (3.5)). Here we display 
the explicit expression for the twist-four DA we have de¬ 
rived: 


<l>( 4 )(x; p 2 ) = —4 J 2 (/x 2 ) ^ x 2 ln(x)(l + 02 + 04) — 5a2 xx 

— 7o 4 xx 


— — -xx + 2 x ln(x) 


1 + -^-xx — -j-(xx ) 2 + (13 — 27x + 18x 2 ) xln(x) 


+ (x —> x) 


( 6 ) 


r 


We estimate within this renormalon-based approach 
error ellipses of the CLEO data and show the results 
of the data processing in Figs. 0 and [ 2 ] As one sees 
from these figures, the general effect, entailed by the in¬ 
clusion of the twist-four contributions by means of the 
renormalon approach, is to shift the error ellipses further 
to the right, somewhat enlarging them as well. One ob¬ 
serves that larger values of 02 and 04 are preferred, while, 
at the considered 1<t error level, the sign of 04 remains 
unchanged and negative. This phenomenon is in agree¬ 
ment with what we had found before in Q], see Table 1 
there, when we varied the twist-four coupling parameter 
S(p 2 ) from smaller to larger values. 

For the new best-fit point (O), we obtain the following 
valupq p beSt fit ( ren )|'..2 \ _ I n 01 „ n j best fit(ren), 2 \ _ 

Vdilltib 0*2 \r^SY/ — 1 U.O-L dilvj. CL ^ Vr^SY/ — 

—0.25. Note that the y 2 value of that best-fit point 
(Xbest fit(ren) = 0-58) is comparable with the one in the 
standard analysis (7,0 (Xbest fit = 0-47). Analyzing the 
structure of the twist-four pion DA, given by Eq. 0 , 


one finds that it has two distinctive peaks and strongly 
enhanced endpoint regions x = 0 , 1 , in accordance to the 
assertions in and |2l|. But, perhaps somewhat sur¬ 
prisingly, including this contribution into the LCSR anal¬ 
ysis of the CLEO data, the characteristics of the best-fit 
pion DA in the endpoint region remain almost the same 
as found before 11 with the use of the asymptotic 
form for <F>^ 4 ) (a;; /u 2 ). Both DAs show endpoint suppres¬ 
sion, much like the BMS model DA extracted from the 
NLC QCD SRs El , though these two approaches are not 
related to each other. Given that the renormalon-based 
twist-four contribution represents an upper bound 23 . 
we can conclude that the endpoint region of the pion DA 
will not receive further enhancement. 

Figure 0 presents a collage of several proposed models 
for the pion DA, listed in Table^J in comparison with the 
CLEO-data constraints in terms of the la ellipses, and 
the recent lattice measurements of Refs. [i3, [lji|. The 
influence of NNLO radiative corrections—only partially 
known 12(1 —can be investigated by means of varying the 
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TABLE I: Estimates at the normalization scale g.Q = 1 GeV 2 for the Gegenbauer coefficients and the reduced inverse moment 
for several model DAs for the pion, based on the first two Gegenbauer coefficients 02 and a a- The designations 
correspond to those used in Fig. [5] Included are also the theoretical constraints derived from LCSRs and NLC QCD SRs, and 
by analyzing [j] the CLEO data llp|. Also shown are the lattice measurements of m and m and the previous transverse- 
lattice result. [Note that the uncertainties on the Gegenbauer coefficients 02 and <24 are correlated. Here, the rectangular limits 
of the fiducial ellipse [V, are shown.] 


DA models/methods 

Symbols 

a 2 (/2o) 

04 (^ 0 ) 


As 

♦ 

0 

0 

0 

ADT [26] 

A 

0.05 

-0.04 

0.01 

BMS best fit point [7] 

+ 

0.31 

-0.35 

-0.04 

Best fit point (ren) 

O 

0.44 

-0.40 

0.04 

PPRWG [27] 

☆ 

0.046 

0.007 

0.05 

SY [25] 

• 

0.27 

-0.22 

0.05 

BMS [12] 

X 

0.20 

-0.14 

0.06 

PR [28] 

V 

0.09 

-0.02 

0.07 

BZ [29] 

▲ 

0.12 

-0.02 

0.10 

Agaev [9] 

0 

0.23 

-0.05 

0.18 

CZ [16] 

■ 

0.56 

0 

0.56 

BF [25] 

□ 

0.44 

0.25 

0.69 

NLC QCD SRs for {£ N ) n [8,J+] 


[0.13,0.25] 

[-0.04, -0.22] 

— 

“Daughter” NLC QCD SR [8,JJ 


— 

— 

0.09 + 0.10 

LCSR analysis of JLab data [3b] 


0.24 + 0.14 + 0.08 

— 

— 

LCSR analysis of CLEO data [^, ^b] 


— 

— 

-0.03 + 0.18 

Transverse lattice [M] 

▼ 

0.08 

0.02 

0.10 

QCDSF/UKQCD lattice [17] 

HI 

0.32 + 0.11 

— 

“0.32 + 0.11” 

Del Debbio et al. (lattice) 

■■ 

0.34 + 0.21 

— 

“0.34 ± 0.21” 


factorization scale, say, in the interval Q 2 /2 < / 1 2 < 2 Q 2 , 
where fj , 2 = = /i^. This has been discussed in jgj and 

it has been shown that it does not alter qualitatively the 
main results. Note that the NNLO corrections within the 
conformal scheme have been quantitatively analyzed in 
|2Cj . Looking more carefully on Fig. [3 one recognizes an 
interesting pattern: with the exception of the BF model 
and the CZ one, all other models steer themselves along 
an “orbit” defined by 02+04 « const. The organizing 
tool behind that pattern is the reduced inverse moment 
(a; -1 )^ of the pion DA, 

( x 1 )w = -—^-1 = a 2 + 04 + ... , (7) 

where the ellipsis stands for still higher Gegenbauer coef¬ 
ficients. This inverse moment plays a special role in the 
description of several form factors of the pion in pertur¬ 
bative QCD. 

The rationale for it is as follows. Though the diversity 
in shape among the models shown in the first column of 
Table [I] in correspondence with Fig. [3 As substantial, 
the associated values for the reduced inverse moment are 
relatively close to each other (see the last column of Ta¬ 
ble 0. Nevertheless, the predictions for the pion-photon 
transition form factor vary quite strongly along this orbit 
as it is evident from Fig. [3 The deep reason for this, is 


that this form factor is not directly proportional to the 
reduced inverse moment (for more details, the interested 
reader is referred to Appendix E of [?] and to 0 )- 

In addition, as it can be shown by using the techniques 
developed in [33|, if the reduced inverse moment is 
in the range of the NLC QCD SRs (see Table [I]), then 
also the hard contribution to the electromagnetic pion 
form factor has the right magnitude—provided one uses 
for all model DAs, the same soft contribution, say, that 
employed in [33, |33j on account of quark-hadron local 
duality. Moreover, the estimate (x~ l )^ M /3 — 1 = 0.24 + 
0.16, obtained in the data analysis of the electromagnetic 
(EM) pion form factor within the framework of a different 
LCSR method in [3(| 0 —though only dependent on 
a 2 —is compatible within errors with the quoted NLC 
QCD SR result. On that basis, one might expect that 
future lattice measurements should yield negative values 
of 04 in order to be in line with the typical value of the 
reduced moment. This would decrease the values shown 
in Table^in quotation marks, recalling that they contain 
only the first Gegenbauer coefficient a 2 . 

We have given a comprehensive compendium of the 
current situation on the constraints imposed on the first 
two non-trivial Gegenbauer coefficients by the analysis of 
the CLEO data 0] in comparison with the most recent 
estimates from modern lattice technology 00. We 
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used in our data analysis the method of light-cone sum 
rules including, on one hand, perturbative corrections up 
to NLO, and, on the other hand, taking into account the 
twist-four contribution to the pion DA, having recourse 
to the renormalon approach. We believe that this inves¬ 
tigation further restricts the available space of values in 
the ( 02 , 04 ) plane for the pion DA, confirming the nega¬ 
tive sign of the Gegenbauer coefficient 04 at the ltr level 
of the CLEO-data analysis without and with the inclu¬ 
sion of the twist-four contribution via the renormalon 
method. The gross profile emerging from this CLEO- 
data analysis is consistent with that we have determined 
independently some years ago from QCD sum rules 
with nonlocal condensates. In fact, all constraints we 
have presented above are consistent with a profile for the 
pion DA that has two humps, but with strongly sup¬ 
pressed endpoints at x = 0 and x = 1. The endpoint- 
enhancement of the renormalon-based twist-four contri¬ 
bution does not change these characteristics. We empha¬ 
size that this type of pion DA leads to a pion-photon 
transition form factor that was found IKI fill to be in 
good agreement with the (low) Q 2 CELLO data (15] as 
well. A full-blown reconstruction of the pion DA in terms 
of its Gegenbauer decomposition may become possible in 
the next years, once the lattice community can provide 
estimates for 04 that seems quite difficult at present. Fur¬ 
ther advances in the calculation of the higher-twist con¬ 
tributions to the pion DA and inclusion of the partial 
NNLO radiative corrections in the standard scheme 01 
will also improve the theoretical precision in assessing 
the pion structure via the CLEO or other data of similar 
precision. In this respect, a combined use of data on the 


pion-photon transition with such on the electromagnetic 
pion form factor, after the upgrade to 12 GeV, of the 
CEBAF machine at JLab, will considerably help to tag 
the pion structure at a still deeper level. 

Two of us (A.P.B. and S.V.M.) are indebted to Prof. 
Klaus Goeke for the warm hospitality at Bochum Uni¬ 
versity, where part of this work was carried out. N.G.S. 



FIG. 2: Collage of several proposed models, compiled in Ta¬ 
ble [I] in comparison with the CLEO data and lattice mea¬ 
surements. The dashed rhombus containing the BZ point 
represents the BZ constraints l2fj . All results are evaluated 
at /x| Y = 5.76 GeV 2 after NLO ERBL evolution. 
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